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Abstract 

In this work we combine the large per-photon optical gradient force with the sensitive feedback of a high 
quality factor whispering-gallery microcavity. The cavity geometry, consisting of a pair of silica disks sep- 
arated by a nanoscale gap, shows extremely strong dynamical backaction, powerful enough to excite giant 
coherent oscillations even under heavily damped conditions (mechanical Q rj 4). In vacuum, the thresh- 
old for regenerative mechanical oscillation is lowered to an optical input power of only 270 nanoWatts, 
or roughly 1000 stored cavity photons, and efficient cooling of the mechanical motion is obtained with a 
temperature compression factor of 13 dB for 4 micro Watt of dropped optical input power. 
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Many precision position measurement devices involve the coupling of mechanical degress of 
freedom to an electromagnetic interferometer or cavity^. Today, cavity-mechanical systems span 
a wide range of geometries and scales, from multi-kilometer long gravitational-wave detectors^ 
to coupled nanomechanical-microwave circuits^. For the sensitive detection and actuation of me- 
chanical motion, each of these systems depend upon "dynamical backaction'^ resulting from 
the position-dependent feedback of electromagnetic wave momentum. Recent work in the optical 
domain has used the scattering radiation pressure force to both excite and dampen oscillations 
of a micro-mechanical resonato r 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! , with the intriguing possibility of self-cooling the 
mechanical system down to its quantum ground-state. 

As has been recently proposecP^-^ and demonstrate d 16 ! 17 i the optical gradient force within 
guided-wave nanostructures can be orders-of-magnitude larger than the scattering force. In this 
work we combine the large per-photon optical gradient force with the sensitive feedback of a 
high quality factor whispering-gallery microcavity. The cavity geometry, consisting of a pair of 
silica disks separated by a nanoscale gap, shows extremely strong dynamical backaction, powerful 
enough to excite giant coherent oscillations even under heavily damped conditions (mechanical 
Q ^ 4). In vacuum, the threshold for regenerative mechanical oscillation is lowered to an optical 
input power of only 270 nano Watts, or roughly 1000 stored cavity photons, and efficient cooling of 
the mechanical motion is obtained with a temperature compression factor of 13 dB for 4 micro Watt 
of dropped optical input power. These properties of the double-disk resonator make it interesting 
for a broad range of applications from sensitive force and mass detection in viscous environments 
such as those found in biology^U^, to quantum cavity-optomechanics in which a versatile, chip- 
scale platform for studying the quantum properties of the system may be envisioned. 

The per photon force exerted on a mechanical object coupled to the optical field within a res- 
onant cavity is given by /zgoM> where goM = d(O c /dx is a coefficient characterizing the disper- 
sive nature of the cavity with respect to mechanical displacement, x. In a Fabry-Perot (Fig. [T^) 
or microtoroid resonator (Fig. [TJ5), the optical force manifests itself as a so-called scattering ra- 
diation pressure due to direct momentum transfer from the reflection of photons at the cavity 
boundar y 120 ! 2 H As the momentum change of a photon per round trip is fixed inside such cavi- 
ties, while the round-trip time increases linearly with the cavity length, the radiation pressure per 
photon scales inversely with the cavity size. In contrast, for the gradient optical force the cavity 
length and the optomechanical coupling can be decoupled, allowing for photon momentum to be 
transfered over a length scale approaching the wavelength of Iigh l 14 ! 15 [ This method was recently 
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employed in a silicon photonic circuit to manipulate a suspended waveguide . However, with- 
out the feedback provided by an optical cavity or interferometer, the optical force only provides a 
static mechanical displacement. 

In the case of a cavity optomechanical system, dynamical backaction can be quantified by 
considering the magnitude of the damping/amplification that an input laser has on the mechanical 
motion. For a fixed absorbed optical input power in the bad-cavity limit (fc £2 m ), the maximum 
rate is given by 



where (o c is the optical cavity resonance frequency, m x is the motional mass of the optomechanical 
system, is the optical power dropped (absorbed) within the cavity, and K = K e / K", is a cavity 
loading parameter (?q, the intrinsic energy loss rate of the optical cavity; K e , the energy coupling 
rate between external laser and internal cavity fields). The effectiveness of the coupling between 
the optical and mechanical degrees of freedom can thus be described by a back-action parameter, 
B = gQ M / (K?tt) c m x ) , which depends upon the motional mass, the per-photon force, and the optical 
cavity g-factor. 

Here we describe the design, fabrication, and characterization of a nano-optomechanical sys- 
tem (NOMS) consisting of a pair of optically thin disks separated by a nanoscale gap. The double- 
disk structure (Fig.[l]:) supports high-<2 whispering-gallery resonances, and provides back-action 
several orders of magnitude larger than in previously demonstrated gradient force optomechani- 
cal systems^E^ (very recent workP^ involving the versatile coupling of external nanomechanical 
elements to the near-field of a high-Finesse microtoroid has realized very strong dynamical back- 
action, although still roughly two-orders of magnitude smaller than in our integrated device). The 
double-disk cavities are formed a silicon wafer on which a multi-layer stack has been deposited 
consisting of silicon dioxide disk layers separated by a sacrificial amorphous silicon (a-Si) layer. 
Formation of the circular disk shape along with supporting fork structures are defined by electron- 
beam lithography and an optimized plasma dry-etch. A dry release of the silica disk layers is 
performed using a highly- selective plasma etch of the Ct-Si intermediate layer and the underlying 
Si substrate (see Methods). The final double-disk structure, shown in Fig. [2j consists of 340-nm- 
thick silica disks separated by a 140 nm air gap extending approximately 6 /im in from the disk 
perimeter (the undercut region). Two different sized cavities are studied here, one large (D = 90 
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FIG. 1: Schematic of the corresponding (a) Fabry-Parot and (b) microtoroid optomechanical cavities, 
(c) Schematic of the double-disk NOMS structure, showing the mechanical flapping motion of the disks. 
FEM-simulated optical mode profiles of the radial component of the electric field for the (d) bonded mode 
at A = 1520 nm and (e) antibonded mode at X = 1297.3 nm. (f) FEM-simulated tuning curve of the bonded 
mode, (g) Optomechanical coupling coefficient and effective length (blue curves) for the bonded mode. 
gOM and Lqm are both well-approximated by exponential functions (red curves). 

/im; Sample I) and one small (D — 54 /im; Sample II) in diameter. The small diameter cav- 
ity structure represents a minimal cavity size, beyond which radiation loss becomes appreciable 

(Qr ~ 10 8 ). 

Finite element method (FEM) simulations of the whispering-gallery optical modes of the 
double-disk structure shows substantial splitting of the cavity modes into even and odd parity 
bonded and anti-bonded modes (Fig.[l£-f). Due to its substantial field intensity within the air gap, 
the bonded mode tunes rapidly with changing gap size as shown in the inset to Fig. [TJ*. The op- 
tomechanical coupling coefficient, g<jM, defined as the derivative of the mode tuning curve, can be 
related to an effective optomechanical coupling length, Lqm through the relation goM = 0) c /Lqm- 
As the two disks are coupled through the evanescent field between them, Lqm decreases expo- 
nentially with disk spacing (Fig. [I]*), reaching a minimum value of 3.8 /im at a resonance optical 
wavelength of X c m 1.5 /im. For the air gap of 138 nm used in this work, the optomechanical 
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coupling is estimated to be goy[/2n = 33 GHz/nm (Lom = 5.8 /im), equivalent to 22 fN/photon. 

The double-disk structure also supports a number of different micro-mechanical resonances, 
ranging from radial breathing modes to whispering-gallery-like vibrations of the disk perimeter. 
The most strongly coupled mechanical resonance is that of the symmetric (i.e., azimuthal mode 
number, m = 0) flapping motion of the disks. Due to the symmetry of the gradient force on the 
two disks, only the 6 /im undercut air gap region is involved in the flapping motion (Fig. [T]: and 
Fig. [2]), with negligible displacement of the interior of the disk. Consequently, the flapping mode 
exhibits an effective motional mass, m x , of only 145 and 264 picrograms for the 54 and 90 /im 
diameter cavities, respectively (see App. [B]). Note that both these values are more than two orders 
of magnitude smaller than commonly used micromirrors and microtoroid g * * ' ' ' ' \ and in 
combination with the large per-photon force, provide a significant enhancement to the dynamic 
back-action parameter which scales as gQ M /m x . 

Optical and mechanical measurements were initially performed at room temperature in a one 
atmosphere nitrogen environment. Fig. [3^ shows the wavelength scan of a large diameter double- 
disk cavity (Sample I). Several radial-order whispering-gallery modes are evident in the spectrum, 
all of them of TE-like polarization and bonded mode character. The fundamental TE-like bonded 
optical mode at X = 1518.57 nm is shown in the Fig. [3^ inset, from which an intrinsic optical 
<2-factor of 1.75 x 10 6 is inferred (see App. [cj. The radio-frequency (RF) power spectrum of the 
optical signal transmitted through the cavity (Fig. |3]5, top panel) exhibits three clear frequency 
components at 8.30, 13.6, and 27.9 MHz corresponding to thermally-actuated resonances of the 
double-disk structure. These values agree well with FEM simulations of the differential flapping 
mode mode (7.95 MHz), and the first (14.2 MHz) and second (28.7 MHz) order radial breathing 
modes (Fig. [3]:). The strong dynamic back-action of the flapping mode (under thermal excitation) 
also produces a broadband spectral background in the RF spectrum with a shoulder at the second 
harmonic frequency (see App. [C]). The correct description of the power spectrum (Fig. |3]5, red 
curve) shows that the flapping mode has a 3-dB linewidth of 2.1 MHz (mechanical Q-factor, 
Qm = 3.95), limited by the squeeze-film process of the nitrogen gas between the disksP. 

Despite the near-unity mechanical quality factor of the flapping mode, the powerful dynamic 
back-action in the double-disk structure provides sufficient compensation of mechanical loss to ex- 
cite regenerative mechanical oscillation. As shown in Fig.[3]l, with an input optical power of 760 
/iW launched at the blue detuned side of the resonance, the induced parametric mechanical insta- 
bility causes the cavity transmission to oscillate over the entire coupling depth with a fundamental 
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FIG. 2: (a) Schematic of the experimental setup for optical testing of the double-disk cavity. The cavity 
input and transmission are both transported through a single-mode silica fiber taper, which is supported 
by two nanoforks for stable operation. A tunable laser source is used to optically probe and actuate the 
double-disk structure, with input power controlled by a variable optical attenuator (VOA) and wavelength 
calibrated by a Mach-Zehnder interferometer (MZI). For experiments performed in a nitrogen environment, 
the cavity transmission is sent directly to the photodetectors, while it is first amplified by an erbium-doped 
fiber amplifier (EDFA) for the experiments performed in vacuum, (b-c) Scanning electron microscope 
images of the 54-jUm double-disk NOMS. False color is used to indicate different relevant regions of the 
device. 



frequency of 13.97 MHz (this value is about 68% larger than the intrinsic mechanical frequency 
due to the optical spring effect 24 ). A zoom-in of the recorded time waveform (Fig. [3^) agrees well 
with our numerical simulation which shows that the gradient force actuates an extremely large (50 
pm) mechanical displacement amplitude, dragging the cavity resonance over more than 10 cavity- 
linewidths and leaving distinctive features of the Lorentzian cavity transfer function. In particular, 
two sequential passes of the cavity resonance across the laser frequency can be seen, along with 
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FIG. 3: (a) Optical transmission spectrum of a large diameter (D = 90 ^m; Sample I) double-disk cavity. 
The inset shows the fundamental TE-like bonded mode at A = 1518.57 nm. (b) Upper panel: optical 
transmission power spectral density (PSD) of a Sample I double-disk in the 1 atm. nitrogen environment 
for P, ' = 5.8 juW. Experimental data in blue, theoretical modeling in red, and detector noise background in 
yellow. Lower panel: transmission PSD of a small diameter (D = 54 /im; Sample II) double-disk cavity in 
vacuum for P = 44 nW. The inset shows a zoomed-in of the spectrum around the fundamental flapping mode 
frequency, (c) FEM simulated mechanical modes indicated in (b). (d) Recorded transmission waveform 
of Sample I for Pi = 0.76 mW. (e) Comparison of experimental (blue curve) and simulated (red curve) 
waveforms, with the corresponding simulated mechanical displacement (green curve). Further details of 
the measurement conditions are given in App. [Dl 



an overshoot and oscillation of the transmitted optical power resulting from the quick release of 
Doppler shifted photons from the cavity. 

The threshold for regenerative oscillation depends sensitively upon the optical input power and 
the average laser-cavity resonance detuning, a map of which can be used to quantify the strength of 
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the dynamic back-action. An estimate of the threshold detuning (A t h), for a given input power, can 
be determined from the abrupt kink in the cavity transmission that marks the onset of regenerative 
oscillation (Fig. [4^ and App. [D]). The detuning dependence of the optomechanical amplification 
coefficient can be lumped into a single detuning function, 

^K A -^)( (A+Q » )2+ © 2 )( (A - fl ™ )2+ ©> (2) 

where K = JQ + K e is the total photon decay rate of the loaded cavity. The right panel of Fig. |4j) 
shows a map of /(A^) versus optical input power for the 90 /im diameter double-disk cavity in the 
heavily damped nitrogen environment. The data in Fig.[4j>, as expected, shows a linear dependence 
of /(A t h) on input power, and is well described in the unresolved sideband regime^ by 

fM-J^rt (3) 
CQ c m x r m K- \i-mj 

where T m = 2.1 MHz is the bare mechanical damping rate of the flapping mode. Fitting of eq. ([3]) 
to the data in Fig.[4j) yields a dynamic back-action parameter of B = 0.061 MHz//iW, correspond- 
ing to an optomechanical coupling factor of ^om/2tt = 33.8 ±0.4 GHz/nm, in good agreement 
with the simulated result of 33 GHz/nm. 

In order to eliminate the squeeze-film damping of the nitrogen environment, measurements 
were also performed in vacuum (P < 5 x 10 -4 Torr). The significantly reduced mechanical 
linewidth in vacuum shows that the flapping mode consists of a small cluster of modes (Fig. [3}), 
bottom panel). As discussed more fully in App. |Dj these modes are a mixture of the lower-lying 
azimuthal modes, coupled together due to deviations in circularity of the undercut region and 
support pedestal. Measurements of the optical spring effect indicates that the optical field renor- 
malizes the cluster of modes, with the lowest-frequency mode at 8.53 MHz transforming into the 
fundamental flapping mode with uniformly distributed displacement along the disk perimeter (the 
rest of the modes decouple from the light field). With an in-vacuum <2m = 4070 (Fig. |3]5, inset), 
the fundamental flapping mode has an extremely low threshold input power for regenerative os- 
cillation. Figure [4^ shows a transmission spectrum when the laser is scanned across the cavity 
resonance. Three different regimes can be clearly seen: (I) transduction of thermal motion, (II) 
onset of optically-driven oscillation, and (III) optically damped motion. The onset of regenera- 
tive oscillation coincides with a frequency shift in the fundamental flapping mode to 9.63 MHz 
as shown in the bottom two panels of Fig. [4^. The left panel of Figure |4j) shows a plot of the 
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FIG. 4: (a) Top panel: Normalized cavity transmission for Sample II in vacuum and Pi = 11 /iW. Blue and 
red traces show the instantaneous and low-pass-filtered signals, respectively. Middle panel: the transduction 
amplitude of the frequency component at 8.53 MHz and its higher-order harmonics. Bottom panel: the 
transduction amplitude of the frequency component at 9.63 MHz and its higher-order harmonics, (b) /(A t h) 
as a function of optical input power. Right Panel: Sample I in a 1 atm. nitrogen environment. Left panel: 
Sample II in vacuum (inset shows the minimum achievable threshold (green arrow)), (c) Spectral intensity 
of the thermally-driven fundamental flapping mode at various input powers, recorded for Sample II in 
vacuum, with a laser detuning of A = — 1.45(fc/2) (inset shows the displacement sensitivity at the highest 
input power with the second optical attenuator removed), and (d) the corresponding effective temperature. 
In (d) the red curve is a fit to the data, the solid green (dashed black) curve is a theoretical curve obtained 
using the estimated B-parameter from the left panel of (b) and the experimental (optimal) detuning of 
A=-1.45(k/2) (a = -(k/2)/V5). 

in- vacuum /(A t h) versus input power, with a measured minimum threshold power of Pi = 267 nW. 
Extrapolation of the experimental data using Eqs. ([2]) and ([3]) to the optimal detuning point shows 
a minimum threshold power of only 40 nW. 

The large mechanical amplification of the double-disk NOMS implies a correspondingly effi- 
cient cooling of mechanical motion on the red-detuned side of the cavity resonance. As shown 



9 



in Fig. |4j: for Sample II in vacuum, the spectral intensity of the fundamental flapping mode de- 
creases dramatically with increased input power, accompanied by a significant broadening of the 
mechanical linewidth. Even for the strongest damping levels, the inset to Fig. Hb shows good sig- 
nal to noise for the transduced motion due to the high displacement sensitivity of the double-disk 
(7 x 10~ 17 m/Hz 1 / 2 , as limited by the background level). A measure of the optical cooling can 
be determined from the integrated area under the displacement spectrurrP^ (see App. [e]). Figure 
|4jl plots the inferred temperaure, T e g, which drops down to 12.5 K for a maximum input power 
of Pi = 11 /iW (Pd = A A /iW). In principle, the effective temperature is related to the optical 
damping rate (r m opt ) through the relation Tq/T^ = 1 +r m /r m op t, where Tq = 300 K is the bath 
temperature. In Fig. [4]d the red curve is a fit of the measured cooling curve using the relation 
To/T e ff = 1 + aPi, whereas the green curve represents the expected cooling curve for the dynamic 
back-action parameter (B = 0.032 MHz//iW) determined from the threshold plot in the right panel 
of Fig. [4}) and the experimental laser-cavity detuning (A = — 1 .45(jf/2)). For comparison, we have 
also plotted (dashed black line) the theoretical cooling curve in the case of optimal laser-cavity de- 
tuning (A = — (k/2)/v / 5). The difference between the two theoretical curves and the measured 
data, along with the limited range of optical input power studied, can largely be attributed to issues 
associated with the limited bandwidth and range of our current cavity locking scheme (a problem 
exacerbated by the very large transduction of even the Brownian motion of the disks). As the 
dashed black curve indicates, technical improvements in the cavity locking position and stability 
should enable temperature compression factors of 20 dB for less than 1 /iW of dropped power. 

The large dynamic back-action of the double-disk cavity, primarily a result of the large per- 
photon force and small motional mass of the structure, opens up several areas of application out- 
side the realm of more conventional ultra-high-<2 cavity geometries. This can be seen by consid- 
ering not only the efficiency of the cooling/amplification process, but also the maximum rate of 
effective cooling/amplification, the scale of which is set by the optical cavity decay rat d 26 * 27 [ In 
the double-disk cavities presented here, the dynamic back-action parameter is B « 0.06 MHz//iW 
for a cavity decay rate of k/2tz « 100 MHz. The combination allows for higher mechanical fre- 
quencies of operation, where the bare damping is expected to scale with frequency, and makes 
possible enormous temperature compression ratios. A quantum mechanical analysis of the optical 
self-cooling process^ 2 ^, indicates that the sideband resolved regime (k < y/32Q. m ) is necessary 
to reduce the phonon occupancy below unity. Having already achieved optical Q-factors in excess 
of 10 6 , and planar silica microdisks having already been demonstrated with Q > lO 7 ^, we expect 
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that further optimization of the double-disk NOMS will be able to extend its operation well into 
the sideband resolved regime. The combination of large dynamic back-action parameter and large 
maximum amplification rate also present intriguing possibilites for sensitive, high temporal reso- 
lution force detection 29 , particularly in heavily damped environments such as fluids for biological 
application ^ 18 ! 19 !. Other application areas enabled by the chip-scale format of these devices include 
tunable photonic s-^EDS optical wavelength conversion^!, and RF-over-optical communication. 
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APPENDIX A: DOUBLE-DISK FABRICATION 

Fabrication of the double-disk whispering-gallery resonator consisted of initial deposition of 
the cavity layers. The two silica disk layers and the sandwiched amorphous silicon (a-Si) layer 
were deposited on a (100) silicon substrate by plasma-enhanced chemical vapor deposition, with 
a thickness of 340 ± 4 nm and 158 ± 3 nm for the silica and a-Si layers, respectively. The wafer 
was then thermally annealed in a nitrogen environment at a temperature of T = 1050 K for 6 
hours to drive out water and hydrogen in the film, improving the optical quality of the material. 
The disk pattern was created using electron beam lithography followed by an optimized C4F8- 
SFg gas chemistry reactive ion etch. Release of the double disk structure was accomplished using 
a SFg chemical plasma etch which selectively (30,000 : 1) attacks the intermediate CC-Si layer 
and the underlying Si substrate, resulting in a uniform undercut region between the disks which 
extends radially inwards 6 /im from the disk perimeter. Simultaneously, the underlying silicon 
support pedestal is formed. The final gap size between the disks was measured to be 138 ± 8 nm 
(shrinkage having occurred during the anneal step). Two nanoforks were also fabricated near the 
double-disk resonator to mechanical stabilize and support the fiber taper during optical coupling; 
the geometry was optimized such that the forks introduce a total insertion loss of only ~8%. 

APPENDIX B: OPTICAL GRADIENT FORCE AND EFFECTIVE MOTIONAL MASS 
1. Optical gradient force in a double-disk NOMS 

As the mode confinement in a double-disk NOMS is primarily provided by the transverse 
boundaries formed by the two disks, the double-disk structure can be well approximated by a 
symmetric double- slab waveguide shown in Fig. [5} For the bonding mode polarized along the e y 
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FIG. 5: Schematic of a symmetric double-slab waveguide, h and xq are the slab thickness and the slab 
spacing, respectively. n c and n s are the refractive indices for the slab and surrounding area, respectively. 



direction, the tangential component of the electric field is given by: 
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where x is the transverse component of the propagation constant inside the slabs and 7 is the field 
decay constant in the surrounding area. They are given by the following expressions: 



where Icq = COq/c is the propagation constant in vacuum and /3 = &on e ff is me longitudinal com- 
ponent of the propagation constant of the bonding mode. n e g is the effective refractive index for 
the guided mode. Accordingly, the tangential component of the magnetic field can be obtained 
through H z = The continuity of E y and H z across the boundaries requires K and 7 to 

satisfy the following equation: 



which reduces to tan Kh = y/x when xq — > 0, as expected. 

The circular geometry of the double disk forms the whispering-gallery mode, in which the 
resonance condition requires the longitudinal component of the propagation constant, /3, to be 
fixed as 2nRf} = 2mK, where R is the mode radius and m is an integer. Thus, any variation on 



(B2) 



Ky[l+tanh(y* /2)] = 



2 - 7 2 tanh(7xo/2)] tan xh, 



(B3) 



14 



the disk spacing x$ transfers to a variation on the resonance frequency (Oq through Eqs. (B2) and 
( |B3[ ), indicating that C0q becomes a function of xq. By using these two equations, we find that the 
optomechanical coupling coefficient, gQM = ^> is given by the general form 



gOM(*o) 



A(n 2 -n 2 ) tan Kh + n 2 x Xsech 2 (^) + 2E, (n 2 yhcsc 2 Kh + 2n 2 ) tan Kh + 



npc 

y 



n 2 c y 



where % = K + 7 tan *h an d £ = 1 + tanh( ^ ) . 

When — > 0, Eq. ( |B4| ) leads to the maximum optomechanical coupling as 



£om(0) 



(Oof 



(B5) 



2j3 2 + 2kln 2 c yh' 

In analogy to Fabry-Parot cavities and microtoroids, the magnitude of the optomechanical coupling 



can be characterized by an effective length, L<jm, defined such that goM = jrj^- Equation (B5) 
infers a minimum effective length 



2 
7 



Y 



A n eff + fc 0KV n eff-^ 



K 



{nln-n 2 f 2 



(B6) 



which is approximately on the order of the optical wavelength Xq. 

Physically, as the two slabs are coupled through the evanescent field between them with am- 



plitude decaying exponentially with slab spacing at a rate y [see Eq. (Bl )], the resulting optome- 



chanical coupling can be well approximated by an exponential function 



gOM(*o) ~£OM(0)e 



-7*o 



(B7) 



where £om(0) is given by Eq. (B5). As indicated by the red curve in Fig. 1(g) in the main text, 



Eq. (B7) provides an excellent approximation for the optomechanical coupling coefficient in a 



double-disk NOMS. Therefore, the approximate effective length, Lom ~ g^ToT ^ ' a S rees we ^ 
with the results simulated by the finite element method [see Fig. 1(g) of the main text], and the 
effective length grows roughly exponentially with the disk spacing. 



2. Effective motional mass for the flapping mode 

With a clamped inner edge and a free outer edge, the mechanical displacement of a double 
disk exhibiting a flapping mode is generally a function of radius (Fig. [6]). What matters the op- 
tomechanical effect, however, is the disk spacing at the place where the whispering-gallery mode 
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whispering-gallery mode 



FIG. 6: Illustration of the disk displacement, xq is the disk spacing in the absence of the optical field, ro 
is the radius of the whispering-gallery mode. x m (ro) corresponds to the effective disk spacing at the mode 
center. r a and rj are the inner and outer radii of the disk region involved in the flapping motion. d(r) is the 
mechanical displacement at radius r. 

is located, as that determines the magnitude of the splitting between the bonding and antibond- 
ing cavity modes. As the mechanical displacement actuated by the gradient force is generally 
small compared with the original disk spacing xq, we can assume it is uniform in the region of 
the whispering-gallery mode and define the effective disk spacing x m (ro) at the mode center, 
where ro is the radius of the whispering-gallery mode. The effective mechanical displacement 
is then given by x e ff = x m (ro) — xq, corresponding to an effective mechanical potential energy of 
E p = m x Qj^xZff/2, where m x is the corresponding effective motional mass and Q, m is the resonance 
frequency of the flapping mode. Note that jc e ff is twice as the real displacement at the mode center 
for a single disk, x e ft = 2d(ro). E p reaches its maximum value when the double disk is at rest at its 
maximum displacement, at which point all of the mechanical energy is stored in the strain energy 
U s . Therefore, E p = U s and the effective motional mass is given by 

2U S _ U s 

mx -al[x m (ro)-xo} 2 ~2aid^r y (B8) 

where both U s and d(ro) can be obtained from the mechanical simulations by the finite element 
method. 

The relationship between the effective mass and the physical mass of the double-disk NOMS 
can be found by examining the mechanical potential energy. With a mechanical displacement d(r) 
for each single disk [Fig. [6), we can find the total mechanical potential energy by integrating over 
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the disk regions involved in the flapping motion: 

E p = r Q. 2 m d 2 (r)C2nrhdr, (B9) 
Jr a 

where C, is the material density, h is the thickness for a single disk, r a and are the inner and 
outer radii of the disk region involved in the flapping motion (see Fig. [6]). Note that E p is the total 
potential energy for the two disks, which is simply two times that of single one because of the 
symmetry between the two disks. As the physical mass of a single disk region involving in the 



flapping motion is given by m p = TiC,h{r\ — r 2 ), using Eq. (B9 ), we find that the effective mass is 
related to the physical mass through the following expression: 

>*»= ,, 2U * a rd 2 {r)dr= * / rd 2 (r)dr. (BIO) 

As the whispering-gallery mode is generally located close to the disk edge (i.e., the mode radius 
tq = 44 jim in a double disk with r& = 45 /im), d 2 (r)/d 2 (ro) <C 1 for most of the region between 



r a and rj,, and Eq. (BIO) shows that m x <^ m p /2. Therefore, the effective mass is significantly less 
than half the physical mass of a single disk region. In practice, the effective mass is much smaller 
than this value because of the real displacement function d(r). For the 90-/im device used in our 
experiment, the effective mass is 0.264 nanogram, only about one fifth of the physical mass of a 
single disk region m p = 1.18 nanogram. The effective mass decreases to 0.145 nanogram for the 
54-/im device, due to the decrease in the disk radius. 



APPENDIX C: LINEAR TRANSMISSION OF AN OPTOMECHANICAL CAVITY 

Unlike other microcavities in which the linear transmission is determined only by the cavity 
loss and dispersion, for the double-disk NOMS, even the small thermal Brownian motions of 
the flapping mode introduce significant perturbations to the cavity resonance due to the gigantic 
opto-mechanical coupling, leading to considerably broadened cavity transmission. Figure [Ja) 
shows an example of the cavity transmission of Sample I. With a small input power of 5.8 fiW 
well below the oscillation threshold, the cavity transmission exhibits intense fluctuations when the 
laser frequency is scanned across the cavity resonance. As a result, the averaged spectrum of the 
cavity transmission (red curve) is significantly broader than the real cavity resonance. A correct 
description of the cavity transmission requires an appropriate inclusion of the optomechanical 
effect, which is developed in the following. 
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FIG. 7: (a) The cavity transmission of Sample I in a nitrogen environment, when the laser is scanned across 
the cavity resonance at 1518.57 nm with an input power of 5.8 |iW. The blue curve is the instantaneous 
signal collected by the high-speed detector and the red curve is the average signal collected by the slow 
reference detector 2. The slight asymmetry in the transmission spectrum is due to the static component 
of mechanical actuation when the laser is scanned from blue to red. The dashed line indicates the laser 
frequency detuning used to record the power spectral density shown in the top panel of Fig. 3(b) in the main 
text, (b) Linear scan of the averaged cavity transmission of Sample I at an input power of 2.9 |iW. The 
inset shows a detailed scan for the bonding mode at 15 18.57nm, with the experimental data in blue and the 
theoretical fitting in red. 

When the optical power is well below the oscillation threshold and the flapping mode of the 
double disk is dominantly driven by thermal fluctuations, the mechanical motion can be described 
by the following equation: 

d 2 x „ dx o Frit) „ 
+ T +a 2 x= J±A : (CI) 

dt z dt m x 
where Q, m , T m , and m x are the resonance frequency, damping constant, and effective mass of the 
flapping mode, respectively. Fj is the Langevin force driving the mechanical Brownian motion, a 
Markovin process with the following correlation function: 

(F T (t)F T (t + T)) = 2m x r m k B T8(T), (C2) 



where T is the temperature and kg is Boltzmann's constant. It can be shown easily from Eqs. (CI ) 



and (C2 ) that the Brownian motion of the flapping mode is also a Markovin process with a spectral 
correlation given by (x(Q.i)x* (Q.2)) = 2nS x (Q.i)8(£li — Q2), where x(Q.) is the Fourier transform 
of the mechanical displacement x(t) defined as x(CL) = f_™x(t)e l£lt dt, and S X (Q.) is the spectral 



18 



intensity for the thermal mechanical displacement with the following form: 

2T m ksT /m x 

Sx{a} = (n2-n2)2 + (nr B )2- (C3) 

The time correlation of the mechanical displacement is thus given by 

(x(t)x(t + x)) = — / S x (Q.)e-' aT dT= (x 2 )p(x) « (x 2 )e- T ^ 2 cosa m x, (C4) 

111 J-oo 

where (x 2 ) = ksT /(m^Q^) is the variance of the thermal mechanical displacement and p(r) is the 
normalized autocorrelation function for the mechanical displacement. 

To be general, we consider a doublet resonance in which two optical fields, one forward and the 
other backward propagating, circulate inside the microcavity and couple via Rayleigh scattering 
from the surface roughness. The optical fields inside the cavity satisfy the following equations: 

(z'Ao - k/2 - igoMx)af + ir\a b + {■s/K e A im (C5) 

(iAo - k/2 - igoMx)a b + ir\a f , (C6) 



dcif 

dt 
da b 
dt 



where cif and a b are the forward and backward whispering-gallery modes (WGMs), normalized 
such that Uj = \aj\ 2 (j = f,b) represents the mode energy. A[ n is the input optical wave, normal- 
ized such that Pi n = \Aj n \ 2 represents the input power. K is the photon decay rate for the loaded 
cavity, and K e is the photon escape rate associated with the external coupling. Ao = CO — (Oq is 
the frequency detuning from the input wave to the cavity resonance and rj is the mode coupling 
coefficient. In the case of a continuous-wave input, Eqs. ( |C5[ ) and ( |C6[ ) provide a formal solution 
of the forward WGM: 

p -\-oo 

af(t) = iJKeA in cos(riT)f(T)e-' 80Ml » x{t - T ' )dr 'dT, (C7) 
Jo 

where f(z) = g(* A o"~ K / 2 ) T represents the cavity response. Using Eq. (C4), we find that the statisti- 
cally averaged intracavity field is given as: 

(a f (t)} = i^K- e A in J^ cos(T]T)/(T) e -§^W, (C8) 
where £ = ^qm^ 2 ) an( ^ M T ) ^ s defined as 

h(x) = jj%(x l -x 2 )dx l dx 2 . (C9) 
Similarly, we can find the averaged energy for the forward WGM as: 

n p -\-oo 

(U f (t)} = KePinJJ o f(x l )f(x 2 )cos(r,x l )cos(r,x 2 )e- £ 2 h ^-^dx l dx 2 

= ^V^l re-^[f c (x)+f:(x)]dx + c.c, (CIO) 
2K K — 2iX\ Jo 
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where fj(z) = e(' A /'~ K '/ 2 ) T (j = c,s), with A c = Ao + rj and A s = Ao — 7J. c.c. denotes complex 
conjugate. 

As the transmitted power from the double disk is given by 



MO = ^ + K e U f (t)+iy/K~ e [A? a/(0 -A in a}(t)] 
the averaged cavity transmission, (T) = (Pr)/Pm, tnus takes the form 



(Cll) 



(r) = i 



K e Kj 



1 



lT]K e 



Ki{K-2iT\) 



e -^>[f c (x) + f:(T)}dT + c.c. 



(C12) 



In the case of a singlet resonance, r\ = and Eq. (C12) reduces to the simple form expression 

(r) = l- — / e- £ 2 h ^[f(T)+f(x)}dx. (C13) 

K Jo 



In the absence of opto-mechanical coupling, goM = and Eq. ( |C13[ ) reduces to the conventional 
form of 



T = 1 



(C14) 



A2 + (K/2)2' 

as expected. 

Using the theory developed above and fitting the experimental averaged cavity transmission 
spectrum, we obtain the optical Q factor of the resonance, as shown in Fig. [7](b) for Sample I. The 
same approach is used to describe the cavity transmission of Sample II, given in Fig. 3(a) of the 
main text. 



1. Power spectral density of the cavity transmission 

Here we provide the derivations of the power spectral density of the cavity transmission in the 
presence of mechanical Brownian motion. We present two theories, one for the linear-perturbation 
regime when the optomechanical effect is small, the other a non-perturbation theory accurate for 
arbitrarily strong optomechanical effect. 

a. The linear-perturbation theory 



If the induced optomechanical perturbations are small, Eq. \C1) can be approximated as 

a fit) « iy/KeA; 



i,„ i COs(7]T)/(t) 
10 



1-i'gOM / x{t-x')dx- 
Jo 



dx. 



(C15) 
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In this case, the transmitted optical field can be written as Aj (t) = Ai„ + iy/K^af(t) « Aq + 5A(t) , 
where Aq is the transmitted field in the absence of the optomechanical effect and 8A is the induced 
perturbation. They take the following forms: 



l-Ke C0S(T]T)/(T)JT 



— A-i n Ao, 



(C16) 



SA(t) = igoMK e A in Jtcos(t]t)/(t) / x(t - x')d%' . (C17) 
Jo Jo 

The transmitted power then becomes P(t) = \A T (t)\ 2 m \A \ 2 +Aq<5A(V) +A 8A*(t). It is easy to 
show that (8A(t)) = and (Pr{t)) = \Aq\ 2 . As a result, the power fluctuations, 8P(t) = P T (t) - 
(Pr(t)), become 

SP(t) « g OM Pin / dxu(x) / x{t - z')dx', (C18) 
Jo Jo 

where u{x) = j'jqcos(7]t)[Aq/(t) -A /*(t)]. By using Eq. (C4>, we find the autocorrelation 
function for the power fluctuation to be 



(8P(t)8P(t + t )} tt£Pf n Jj^ dTidT 2 u(Tl)u(T 2 )\]/(tQ,T h T2) 

where ty/(?o, Ti, T2) is defined as 

V/(? ,Ti,t 2 ) = / dx[ dt 2 p(to + t[-T 2 ). 
Jo Jo 



(C19) 



(C20) 



Taking the Fourier transform of Eq. (C19), we obtain the power spectral density Sp(Q.) of the 
cavity transmission to be 



S P (Q) « g 2 0M PlH(Cl)S x (Q) 



(C21) 



where S x (Cl) is the spectral intensity of the mechanical displacement given in Eq. (C3) and H(Q.) 
is the cavity transfer function given by 

2 



H(Q.) 



1 

ay 



u(x)(e i£lT -l)dx 



(C22) 



In the case of a singlet resonance, the cavity transfer function takes the form: 



H(Q) 



(C23) 



[A 2 + (k/2) 2 ] 2 [(Ao + ^ 2 + (k/2)2] [(Aq - a) 2 + (k/2) 2 } ■ 

In most cases, the photon decay rate inside the cavity is much larger than the mechanical damp- 
ing rate, K 3> T m . For a specific mechanical mode at the frequency Q. m , the cavity transfer function 
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can be well approximated by H(£l) H(Q. m ). In particular, in the sideband-unresolved regime, 
the cavity transfer function is given by a simple form of 



H 



[A2 + (K/2)2] 



4 • 



(C24) 



Therefore, Eq. (C21 ) shows clearly that, if the optomechanical effect is small, the power spectral 
density of the cavity transmission is directly proportional to the spectral intensity of the mechanical 
displacement. 



b. The non-perturbation theory 



The situation becomes quite complicated when the optomechanical effects are large. From 



Eq. (Cll), the autocorrelation function for the power fluctuation of the cavity transmission, 
8P(t) = P T (t) - (Pt), is given by 

(8P(h)8P(t 2 )) = ^(U fl U f2 )-lQ((Ala fl -A in a} 1 )(Ala f 2-A in a} 2 )) 

+ ikJ 1 [(Ufi (A* n a f2 -A in a* f2 )) + (Uf2 {A* n a f i -A in a* fl )}] 

- [K e (U f } + iy/^ e (A* n (a f }-A m (a}})] 2 , (C25) 



where Ufj = Uf(tj) and afj = cif(tj) (j = 1,2). Equation (C25) shows that the autocorrelation 
function involves various correlations between the intracavity energy and field, all of which can 



be found using Eqs. (C4) and (C7). For example, we can find the following correlation for the 
intracavity field: 



({ A *n a fi - A ina}i) {Aj n a f 2-A in a* f2 )} 

dT^CiCie-^+^ime-^ + Me^ + c-c.]. 



■ K ePfn 



(C26) 



where, in the integrand, Cj = cos(r/Tj), hj = h(tj), fj = f(Xj) (with j = 1,2), and y = y/fa 



?i , Ti , T2) - h{x) and yfe — h, Ti, %%) are given by Eqs. (C9) and (C20), respectively. 



Equations ( |C9[ ) and (C20) show that h{x) and yf{t 2 —hi^li^l) var y with ti me on ti me scales 
of \/Q. m and 1 /Y m . However, in the sideband-unresolved regime, K^$>T m and K Q. m . As the 
cavity response function /(t) decays exponentially with time at a rate of ?c/2, the integrand in 



Eq. (C26) becomes negligible when X\ 3> 2/k or T 2 2/k. Therefore, y{t 2 — t\,T\,T 2 ) can be 
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well approximated as 



1 



27t(x 2 } 

X\X 2 

2n(x 2 ) J- 



/-|-oo 
S x (Q.)e- ia ^- h Ua = xix 2 p{t 2 -h). (C27) 
-00 



Similarly, h(x) ~ T , since /z(t) = ty/(0, T, t). Therefore, Eq. (C26) becomes 



(Mz«/l -Mna* n ) (Ala f2 -A in a* f2 )) « -K e if„4>(Af,CiC 2 ) : 



(C28) 



where A? = f 2 — *i and <I>(A?,CiC 2 ) is defined as 



^.C,C 2 )= /jf JT 1 JT 2 C 1 C 2 e-§(^ 2 )[ /l/2 e-^^+/ 1 /|e eT '^P +c . c .] 5 (C29) 



with p = p ( A? ) . Following a similar approach, we can find the other correlation terms in Eq. ( C25 ) . 



Using these terms in Eq. ( |C25 ), we find that the autocorrelation function of the power fluctuations 
is given by 

(SP(ti)8P(t 2 )) « K$PfM&,Oi0 2 ) - [Ke(U f ) +i^/Ke (A? <«/) -A in (a}})] 2 , (C30) 



where Oj = o{Xj) (j = 1,2) and <t(t) is defined as 

K e (K 2 + 2r] 2 



o(x) = 



1 



k(k 2 + 4t] : 



cos(7]t) + — j — ^-sin(77T). 
v ' ; K 2 + 4rj 2 v ' ; 



(C31) 



Moreover, Eq. (C8 ) and (CIO) show that, in the sideband-unresolved regime, (a/) and (Uf) are 
well approximated by 



(fl/(0) ~ iV*eAin / COs(TJT)/(T)e 2 T <*T, 

Jo 

/• /»-|-oo r 

(Uf(t)) « K^jf/ /(T 1 )r(T 2 )cOs(T]T 1 )cOs(T 7 T 2 ) e -§( T '- T2 )' 



dx\dx 2 . 



(C32) 
(C33) 



Therefore, we obtain the final term in Eq. (C30) as 

^([/ / ) + i v ^(4(fl / )-A«„(4))^-)cA //~a(T)[/(T)+r(T)] e -§ T2 dT. 







(C34) 



Using this term in Eq. (C30), we obtain the final form for the autocorrelation of the power fluctu- 
ations: 



(8P{h)8P{t 2 )) « K^<J>(A^CTiCT 2 )-4>(oo,CT l0 2)]. 



(C35) 
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It can be further simplified if we notice that the exponential function e ±eT i T 2p(A<) m gq (£29 ) can 
be expanded in a Taylor series as 



,±£TiT 2 p(A/) 



»=o 



(C36) 



Substituting this expression into Eq. (C29) and using it in Eq. (C35 ), we obtain the autocorrelation 
function for the power fluctuation in the following form 



(8P(t)8P(t + t ))^^Pl^ 



£ n p n (t ) 



n=l 
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g:+(-i)"g„ 



(C37) 



where G„ is defined as 



G n = / r n o(r)f(T)e-^ dr. (C38) 
Jo 

In the case of a singlet resonance, rj = and o(x) simplifies considerably to a = Jq/jc The 
autocorrelation function for the power fluctuation is still described by Eq. (C37 ). 

In general, the power spectral density of the cavity transmission is given by the Fourier trans- 



form of Eq. (C37): 



S P (Q) = K 2 e Pl £ 



£ n S n (Q.) 



g:+(-i)' i g„i 



where S n (Ci) is defined as 



p'\x)e ia% dT. 

-co 



(C39) 



(C40) 



Eq. (C3 ) shows that the spectral intensity of the mechanical displacement can be approximated by 
a Lorentzian function, resulting in an approximated p(r) given as p(r) m e m ' T >> cos 



Eq. ([C4|>]. As a result, Eq. (|C40|) becomes 



SJO) 



nT„ 



1 11 

2P J kl(n - k) ! (nr m /2) 2 + [(2k - n)Q. m + £2] 2 ' 



(C41) 



Combining Eq. (C39 ) and (C41 ), we can see that, if the optomechanical coupling is significant, the 
thermal mechanical motion creates spectral components around the harmonics of the mechanical 
frequency with broader linewidths. As shown clearly in Fig. 3(b) of the main text, the second 
harmonic is clearly visible. In particular, if the fundamental mechanical linewidth is broad, various 
frequency components on the power spectrum would smear out, producing a broadband spectral 
background, as shown in the top panel of Fig. 3(b) in the main text for Sample I. This phenomenon 
is similar to the random-field-induced spectral broadening in nuclear magnetic resonance^ and 
atomic resonance fluorescence 32 . 
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The theory developed in this section can be extended easily for the case with multiple mechan- 
ical frequencies. In this case, power spectrum only only exhibits harmonics of each mechanical 
frequency, but also their frequency sums and differences. As shown in the bottom panel of Fig. 3(b) 
of the main text, the frequency components near MHz is the differential frequencies and those 
near 18-20 MHz are the second harmonics and sum frequencies. 



APPENDIX D: REGENERATIVE OSCILLATION AND MECHANICAL RESONANCES 



1. Numerical simulation of optomechanical oscillations 



The optomechanical oscillations are simulated through the following coupled equations gov- 
erning the intracavity optical field and mechanical motions, respectively: 

da K 

— = (zA - — - lgOMX)a + 1^/KeAin, (Dl) 
d 2 x „ dx n Frit) FJt) 

+ r m — + afe = + (D2) 

dt L dt m x m x 

where we have counted in both the thermal Langevin force Ft and the optical gradient force 
F = — for actuating mechanical motions. 



2. Mapping the threshold detuning 

Figure [8] shows an example of the cavity transmission of Sample I. The mechanical flapping 
mode starts to oscillate when the input laser frequency is scanned across a certain detuning. Within 
this detuning value, the same magnitude of optomechanical oscillation is excited over a broad 
range of laser blue detuning. The intense transmission oscillations cover the entire coupling depth, 
leaving an abrupt kink on the transmission spectrum. The coupling depth at the kink point, 
corresponds to the threshold coupling at the given power level, from which we can obtain the 
threshold frequency detuning A^. 



3. Flapping modes with various azimuthal mode numbers 

Because of the extremely short round-trip time of the cavity mode, the optical wave is sensitive 
only to the variations of averaged disk spacing around the whole disk. As a result, the optome- 
chanical coupling for the fundamental flapping mode, which has flapping amplitude uniformly 
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FIG. 8: Scan of the cavity transmission of Sample I at an input power of 0.76 mW, with the instantaneous 
and averaged signals shown in blue and red, respectively. The dashed line indicated the laser frequency 
detuning used to record the time-dependent cavity transmission given in Fig. 3(d) in the main text. 

distributed around the disk perimeter, is maximum, but is nearly zero for flapping mode with 
high-order azimuthal mode numbers. However, due to the asymmetry in practical devices, the 
net variations of averaged disk spacing induced by the high-order flapping modes (with azimuthal 
mode number > 1) is not zero, and their thermal motion is visible in the transmission power spec- 
trum. In general, their optomechanical coupling is weak and does not provide efficient dynamic 
back action. 



APPENDIX E: COOLING OF THERMAL MECHANICAL MOTION 



1. Spectral intensity of optically damped thermal mechanical motion 



In general, the optomechanical effect is governed by Eqs. (Dl ) and (D2 ). However, the optome- 
chanical effect during mechanical cooling is well described by linear perturbation theory since the 
thermal mechanical motions are significantly suppressed. The intracavity field can thus be ap- 
proximated as a(t) ~ ao(t) + 8a(t), where an is the cavity field in the absence of optomechanical 



coupling and da is the perturbation induced by the thermal mechanical motions. From Eq. (Dl ), 
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they are found to satisfy the following equations: 

dao 

dt 
dda 



(iA -K/2)ao + iy/K~ e A in , (El) 



?Aq — K/2)8a — igoMxao- (E2) 



dt 

In the case of a continuous-wave input, Eq. dETl) gives a steady-state value given as: 



«o = -t^ — 7T-, (E3) 

K/2 — iAq 

and Eq. ( |E2| ) provides the spectral response for the perturbed field amplitude, 

V ; i(Ao + a)-K/2' 

where 8d(Q.) is the Fourier transform of 8a(t) defined as 8a(Q) = /+J 8a{t)e iSlt dt. Similarly, 
3c(n) is the Fourier transform of jc(?). 

The optical gradient force, i 7 ,, = — , is given by 

F (t) = r|ao| 2 + a^a(0+a 5a*(0] . (E5) 

The first term is a static term which only affects the equilibrium position of the mechanical motion, 
and can be removed simply by shifting the zero-point of the mechanical displacement to the new 
equilibrium position. Therefore, we neglect this term in the following discussion. The second and 



third terms provide the dynamic optomechanical coupling. From Eq. (E4), the gradient force is 
given by the following equation in the frequency domain: 

p rry) - 2 goi>ol 2A o^) Ag - fl 2 + (k/2) 2 + iKQ 

(A ] coo [(A + ^) 2 + (k/2) 2 ][(A -^) 2 + (k/2) 2 ]- 1 ; 

As expected, the gradient force is linearly proportional to the thermal mechanical displacement. 



Equation ( |D2[ ) can be solved easily in the frequency domain, which becomes 

(Q 2 - a 2 - iT m Q)x = — + — . (E7) 



Equation (E7) together with (E6) provides the simple form for the thermal mechanical displace- 
ment, 

^ = S(ow-tf-ffyj' (E8) 
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where QL and Y' are defined as 



(o! m ) 2 = ai + 



2 ^oMl a o| 2 A 



A 2 - ft 2 + (k/2)' 



m x (Qo [(Aq + Q)2 + (k/2) 2 ] [(Aq - O) 2 + (k/2) 2 ] 



2 So M l a o| 2 Ao 



A 2 -^ 2 +(k/2)^ 



m x coo [(Aq + Q m ) 2 + (k/2) 2 ] [(Aq - On) 2 + (k/2) 2 ] 



(E9) 



r 



2 ^OMl fl o| 2 KAo 



m x coo [(Aq + Q) 2 + (k/2) 2 ] [(Aq - a) 2 + (k/2) 2 ] 



2 <?omM 2 kA 



m x coo [(A + C2 m ) 2 + (k/2) 2 ] [(Aq - n m ) 2 + (k/2); 



(ElO) 



Equations (E8)-(E10) show clearly that the primary effect of the optical gradient force on the 
mechanical motion is primarily to change its mechanical frequency (the so-called optical spring 



effect) and energy decay rate to the new values given by Eqs. (E9) and (ElO). The efficiency of 
optomechanical control is determined by the figure of merit gQ M /m x . On the red detuned side, the 
optical wave damps the thermal mechanical motion and thus increases the energy decay rate. At 
the same time, the mechanical frequency is modified, decreasing with increased cavity energy in 
the sideband-unresolved regime. 



Using Eqs. ( |C2| ) and (jE8j), we find that the spectral intensity of the thermal displacement is 

2T m k B T/m x 



given by a form similar to Eq. (|C3|): 



S X (Q) 



which has a maximum value S x (Qf m 



m x (Cl' m r m f 

placement is equal to the area under the spectrum, 



[(^) 2 -a 2 F+(^) 2 ' (E11) 

2T m k B T rpj i£ var j ance Q f thermal mechanical dis- 



((Sx) 2 ) 



1 

In J- 



S x (Q)dQ 



k B TT n 



(E12) 



Cooling the mechanical motion reduces the spectral magnitude and the variance of thermal dis- 
placement. 



2. Effective temperature of the cooled mechanical mode 

For a mechanical mode in thermal equilibrium, the effective temperature can be inferred from 
the thermal mechanical energy using the equipartition theorem: 

k B T ef f = m x (n'J 2 ((8x) 2 ). (E13) 
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The area under the displacement spectrum thus provides an accurate measure of the effective 
temperature. In practice, fluctuations on the laser frequency detuning may cause the mechanical 
frequency and damping rate to fluctuate over a certain small range [Eq. (E9) and ( |E10[ )], with 



a probability density function of p(Q! m ). As a result, the experimentally recorded displacement 
spectrum is given by the averaged spectrum 

S X (Q) = J S x (Q)p(Q' m )dQ' m , (E14) 



where S x (£l) is given by Eq. (Ell ) and we have assumed / p(Q.' m )d£l' m = 1. The experimentally 
measured spectral area is thus 

— J ^ s x (a)da = J {{8xf) P {Q! m )dQ! m = ((dx) 2 ). (Ei5) 

Therefore, the integrated spectral area obtained from the experimental spectrum is the averaged 
variance of thermal mechanical displacement, from which, according to the equipartition theorem, 
we obtain the effective average temperature 

k B T eff = m x (U' m ) 2 ((5x) 2 ), (E16) 

where Q.' m = J Q! m p{Q! m )dQ! m is the center frequency of the measured displacement spectrum 
S x (Cl). Compared with the room temperature, the effective temperature is thus given by 



r eff (cQ 2 ((8x)*) 
To a 2 m ((dx)\ 

where ((8x) 2 ) is the displacement variance at room temperature, given by the spectral area at Tq. 
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